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Abstract 

The efficiency of an encounter-controlled reaction between two independently-mobile reactants 
on a lattice is characterized by the mean number (n) of steps to reaction. The two reactants 
are distinguished by their mass with the "light" walker performing a jump to a nearest-neighbor 
site in each time step, while the "heavy" walker hops only with a probability p; we associate p 
with the "temperature" of the system. To account for geometric exclusion effects in the reactive 
event, two reaction channels are specified for the walkers; irreversible reaction occurs either in a 
nearest-neighbor collision, or when the two reactants attempt to occupy the same site. Lattices 
subject to periodic and to confining boundary conditions are considered. For periodic lattices, 
depending on the initial state, the reaction time either falls off monotonically with p or displays 
a local minimum with respect to p; occurrence of the latter signals a regime where the efficiency 
of the reaction effectively decreases with increasing temperature. Such behavior can also occur 
when one averages over all initial conditions, but can disappear if the jump probability of the light 
walker falls below a characteristic threshold value. Even more robust behavior can occur on lattices 
subject to confining boundary conditions. Depending on the initial conditions, the reaction time 
as a function of p may increase monotonically, decrease monotonically, display a single maximum 
or even a maximum and minimum; in the latter case, one can identify distinct regimes where 
the above-noted inversion in reaction efficiency can occur. We document both numerically and 
theoretically that these inversion regions are a consequence of a strictly classical interplay between 
excluded volume effects implicit in the specification of the two reaction channels, and the system's 
dimensionality and spatial extent. Our results highlight situations where the description of an 
encounter-controlled reactive event cannot be described by a single, effective diffusion coefficient. 
We also distinguish between the inversion region identified here and the Marcus inverted region 
which arises in electron transfer reactions. 
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1. INTRODUCTION 



Understanding the interplay between diffusion-controlled dynamics and the geometrical 
constraints imposed by the substrate or topology of the host medium is a problem of seminal 
importance today. For example, with recent advances in nanotechnology, one can now 
visualize and manipulate individual atoms and molecules, and design specific geometries 
on atomic length scales |ll Q, The geometrical characteristics of the latter (i.e., size, 
dimensionality and boundary conditions) are found to influence the efficiency of diffusion- 
controlled reactive processes and trapping events governed by specific interactions. 

After the seminal papers by Ovchinnikov and Zeldovich |3| and Toussaint and Wilczek 
documenting an anomalous decrease in the efficiency of the two-species annihilation reaction 
in low dimensions due to reactant segregation effects, other types of encounter-controlled 
processes between unlike species have also received much attention in the literature. In 
particular, the so-called trapping problem and target problem have been extensively stud- 
ied j^, 0]. Such processes are ubiquitous in nature and the underlying theory finds wide 
application in domains such as trapping of charge carriers [1], exciton annihilation and spin 
relaxation processes quenching of delocalized excitations, recombination and coagula- 
tion processes [l^, predator-prey models [^^[iH etc. In the simplest version of the trapping 
problem, diffusing A particles are annihilated upon encounter with immobile traps T fol- 
lowing the scheme A + T —>■ T. The target problem may be thought of as dual to the 
trapping problem jl^.[l3|. i.e., here randomly moving A particles are the ones which destroy 
immobile targets T upon encounter according to the reaction A + T ^ A. The case where 
both A particles and traps diffuse with different diffusivities i s sig nificantly harder to deal 
with and few analytic solutions are known 0, 0, 0, 0, |l^ Il9| ^ With a few exceptions, 
most of the previous work on such variants of the trapping and the target problem has been 
focused on the computation of the survival probability S(t) for the annihilated species and 
the long-time asymptotics of this property as a measure of the reaction's efficiency. 

On the other hand, a number of papers 2^, 21, i^j have taken up the problem of char- 



acterizing the efficiency of two-particle, encounter-controlled reactions on a lattice. The 
type of reaction considered in those references can also be interpreted as a space- and time- 
discrete version of the trapping problem or the target problem as defined above. In contrast 
to previous work, the quantity used to gauge the efficiency of the process was the initial- 
condition-averaged mean number of time steps to reaction (n), whose inverse is a measure 
of the time scale of the reaction and hence the efficiency of the reaction-diffusion event; 
in particular, in the large lattice limit the inverse converges (in magnitude) to the small- 
est eigenvalue of the associated master equation ^. Basically two types of situations were 
considered: a) an asynchronous scenario in which a "light" walker (also termed "walker 1" 
hereafter) performs nearest-neighbor jumps at each time step while a "heavy" walker (walker 
2) remains stationary at a given site, thereby playing the role of an immobile trap and b) 
a situation where both walkers synchronously perform symmetric nearest-neigbor jumps on 
the lattice. Synchronous dynamics turns out to be more efficient in general, but significant 
exceptions are found for sufficiently small lattices. 

If one assumes, for instance in the framework of Arrhenius theory, that the mobility of the 
heavy walker grows with increasing system temperature, one realizes that such exceptions 



^ The counterpart (t) of (n) in the diffusive Umit is connected with the survival probabihty via the relation 
{t) = S{t)dt, (see e.g. ^6], pages 173-174) 
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are quite remarkable. In classical diffusion theory, as formulated by von Smoluchowski and 
Einstein, the diffusion coefficient is a monotonically increasing function of the temperature. 
Thus, in a continuum approximation where one identifies an effective diffusion coefficient the 
efficiency of an encounter-controlled reaction should increase with increasing temperature. 
Deviations from this expectation emphasize the fact that some properties (particularly for 
small systems) are not correctly described by an approach valid for large systems in the 
thermodynamic limit. It is precisely these properties that are of special importance in 
applications of nanotechnology. 

To understand in a more detailed way the role of temperature in influencing the effi- 
ciency of diffusion-controlled events in small systems, the approach taken in |2ol . 21, 2^ 



was generalized in [23[ as follows: a jump probability p of walker was defined such that 
p = corresponded to purely asynchonous behavior and p = 1 corresponded to purely 
synchronous behavior of the two diffusing particles ^. Interpreting the parameter p as a 
measure of the "system temperature", one would intuitively anticipate that an increase in 
temperature should lead to an increase in the efficiency of reaction. For the model studied in 
[is! and this paper, this expectation was confirmed for certain initial conditions. However, 
already in ID, it was shown that for properly chosen initial conditions the model leads to 
the onset of a parameter region where (n) increases monotonically with p, implying that the 
efficiency of the reaction decreases with increasing temperature. 

In what follows, we shall term this region "inverted" by analogy with the well-known 
decrease of the reaction rate with increasing temperature predicted by Marcus for electron 
transfer reactions (j^ . lisl 2^ in certain regimes of parameter space (reaction coordinate, 



iter space (re 



reorganization energy) and later confirmed experimentally [23, |28|, |29| . However, we em- 
phasize that our use of this terminology is purely symbolic, since in the cases with which 
we shall be dealing here the anomalous decrease in efficiency is due to a stricly classical 
interplay between the system geometry and an "interaction zone" (to be specified in Sec. ^ 
in binary collisions, rather than a quantum-mechanical analysis of curve crossing between 
reactant and product states (where Marcus showed that there are regions in which the reac- 
tion probability can actually decrease with increase in "driving force" of the reaction). As 
we shall show, when the size of the "interaction zone" of the two reactants becomes of the 
same order of magnitude as the system size, anomalous behavior is found. 

To explore the relevance of our work to possible experimental realizations, it will be 
necessary to study in detail the extent to which our "inversion region" is robust with respect 
to system dimensionality and/or boundary conditions. For example, the difference equation 
approach developed in [23j for ID periodic lattices cannot be extended to treat the case of 
higher dimensions or non-translationally invariant lattices. Accordingly, we apply here the 
theory of finite Markov processes to investigate the previously-noted, more general cases in 
detail. Further, we shall consider cases where the jump probability of the "light walker" 
is restricted to values less than one (and will find that the "inversion region" disappears 
when the jump probability falls below a threshold value). Finally, we shall propose (several) 
rationale to explain the quantitative behavior uncovered in our study. 

The plan of the paper is then the following: Sec. |21 gives a general definition of the model 
and shows how in the framework of Markov chain theory the mean reaction time and its 



^ Note that, even though we have chosen to interpret the difference in the niobihty of both walkers for p < 1 
as a difference in their masses, other interpretations are also legitimate. In nonequilibrium system one 
could e.g. also speak of a difference in their internal energies. 



3 



variance can be extracted from the so-called fundamental matrix (the size of this matrix 
can be greatly reduced by proper use of the system symmetry). In Sec. El we deal with 
the case of a ID lattice, thereby comparing the results for the periodic and the confining 
case. This is done both for specific initial conditions and at a coarse grained level (average 
over a homogeneous set of initial states). Sec. E] treats the case of a 2D lattice, both for 
periodic and confining boundary conditions. Sec. discusses the case where the light walker 
is allowed to jump with probability less than one at each time step. Sec. El gives a physical 
interpretation of the main results. Finally, Sec. [71 summarizes the conclusions and discusses 
possible extensions of the model. The appendices A and B contain tables with detailed 
analytic expressions for the reaction time in the general case treated in Sec. El 

2. GENERAL FRAMEWORK 

2.1. Definition of the model and use of symmetry 

Consider two reactants (random walkers labeled, respectively, 1 and 2) performing a time 
discrete random walk on an iV-site lattice with dimensionality and boundary conditions to 
be specified later. At each time step, walker 1 performs a symmetric nearest-neighbor jump 
with probability one, while walker 2 either performs a symmetric nearest-neighbor jump 
with probability p (synchronous step) or remains at the same site (asynchronous step). Any 
time the walkers jump on the same site or exchange positions by crossing each other, the 
reaction takes place instantaneously, and the process is terminated. Our main goal will be 
to clarify how the p-dependance of the mean number of time steps to reaction changes with 
the boundary conditions and the dimensionality of the system and to ellucidate the role of 
the initial conditions in this context. 

The stochastic motion of the two walkers in the above model can be considered as an 
absorbing Markov chain. In order to characterize the transient states of the chain let us first 
label the lattice sites from 1 to N. In the case of a ID lattice this can e.g. be done from 
left to right (see Fig. left), while in the case of a 2D lattice one can number the lattice 
sites, say, from left to right and from top to bottom (see Fig. HI right). We shall consider 
two types of boundary conditions, i.e. periodic and confining. By confining boundaries (also 
sometimes termed 'reflecting' in the literature) we mean that whenever either of the two 
walkers attempts to 'exit' the lattice from a boundary site under the dynamics prescribed 
above, it is reset to the same position. 

Each transient state can thus be expressed by a pair of positive integers Xi and X2 
representing the instantaneous position of walker 1 and walker 2 in terms of the above 
lattice site coordinates. Let us hereafter denote by and X2 the initial values of these 
coordinates. Thus, there are N{N — 1) possible pairs of "non-reactive" initial two-particle 
configurations (x5,X2). Each of these configurations has an associated local reaction time 
(n)(2,o ,^.0), whereby the angular brackets denote the average over a uniform ensemble of joint 
random walk realizations. In addition to these quantities, one can construct an expression 
for the global reaction time (n), i.e., the average over a uniform set of initial conditions 
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In the present work we shall make an extensive use of symmetry to simplify the cal- 
culation of (n). Let us e.g. consider the case of a ID lattice with = 4. The set of 
initial configurations obtained by assigning different values to and can be represented 
pictorially as follows: 
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, where the "1(2)" represent sites occupied by walker 1(2) and the "0" represent empty 
lattice sites. 

In the periodic case, the state of the system is fully characterized by the distance d = 
min(|xi — X2I, — |xi — X2I) between the walkers. Thus, the above 12 states can be lumped 
into the first 2 states (2,1) and (3,1), respectively numbered as State 1 and State 2 in the 
short-hand notation to be used in the following. 

Let us now consider the case of a confining lattice. In addition to the first two states, 
one has 4 additional symmetry-distinct states, e.g. (4,1), (3,2), (4,2) and (4,3). In our 
short-hand notation we thus have: 



State 1 
State 2 
State 3 
State 4 
State 5 
State 6 



2-1-0-0 
2-0-1-0 
2-0-0-1 
0-2-1-0 
0-2-0-1 
0-0-2-1 



This numbering procedure for the symmetry-distinct states can now straightforwardly be 
extended to ID confining lattices with arbitrary A^, i.e. by first setting X2 = 1 and letting 
Xi vary from 2 to A^, then increasing X2 by one unit and letting xi vary from 3 to A^, and so 
on. One thus obtains A^(A^ — l)/2 states characterized by xi > X2- In particular, the states 
for the cases A^ = 5 and A^ = 6 to be dealt with in 13.21 respectively read 
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In the case of a periodic lattice, it is sufficient to consider tfie first [A^/2] states, where 
[•] is the floor function (largest integer smaller than the argument). In this case, the state 
number in the short-hand notation is identical with the interparticle distance d. 

In order to compute the global reaction time (n) from the local ones, one constructs 



{n) = J2Mn)i, (2) 
1=1 

where Ng is the number of symmetry distinct states (respectively equal to [N/2] and 
N{N — l)/2 for a periodic and for a confining ID lattice), {n)i is the reaction time when 
both walkers start from a given symmetry-distinct initial state i, and Wi are appropiately 
chosen weights. The latter are easily calculated, e.g. for the periodic ID lattice one has 



for odd values of N and 
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W[N/2] = (4) 

for even values of N. These expressions reflect the fact that for even lattices the statistical 
weight of the initial condition with maximum distance [N/2] between the walkers is half 
the weight of the other initial states (which appear twice when taking the average over a 
homogeneous ensemble of initial conditions), while in the odd lattice case all initial states 
have the same weight. 



2.2. Markovian approach 

In the framework of Markov chain theory the mean number of time steps (n)j before 
termination of a random walk starting from state i is given by the elements of the so-called 
fundamental matrix. This matrix reads 

= °° =k = = 

N^J^Q =(I-Q)-\ (5) 

fc=0 

where Q is the matrix whose elements qnm are the transition probabilities between two 
transient (non-reactive) states n and m. The matrices Q and N thus have the dimension 
Ns X Ng. Denoting by nnm the elements of the fundamental matrix, one then has js^l 

{n)i = ^nij. (6) 
i=i 

For small lattices the analytical computation of the transition matrix Q in terms of 
the jump probability p as well as the subsequent derivation of the fundamental matrix is 
straightforward. The local and the global reaction times are obtained from the matrix 
elements via eqs. (jH)) and (0). The variance {v)i = — (^)f can also be extracted from 
the fundamental matrix. One has bC 



{v),= (^(2N-I)T-T,,j^, z = l,...Ns. (7) 
In this equation, the A^s-dimensional vector T is given by 

T = N (8) 

where ^ is the A^g-dimensional vector with all entries equal to 1. The vector T^g is obtained 
from T by squaring all entries. The average (v) over the initial conditions is obtained by 
evaluating the weighted sum 

Ns 



{v) = J2w^{v)i. (9) 



1=1 



7 



3. ID LATTICE 



3.1. Periodic boundary conditions 

The Markov method can be used to recover results reported earher [i^ using an approach 
based on backward difference equations. The first (rather remarkable) result obtained at 
this level is that the local reaction times (n)j (where the subscript i is by construction equal 
to the initial interparticle distance d^) do not always decay monotonically with the mobility 
p, as one might expect from an analysis of the relative motion of the two walkers using a 
continuum approach. Indeed, monotonic decay of the (n)j's only takes place for even values 
of i, whereas for odd values a local minimum is observed (see Fig. |2I). This minimum shifts 
to higher p values with increasing initial interparticle distance = i. The existence of such 
a minimum means that there is an anomalous inverted region in which > 0, implying 
that the efficiency of the reaction decreases with increasing temperature. Even though this 
behavior holds both for even and odd values of the minima are rather shallow in the odd 
lattice case and are therefore lost when one takes the average over the initial conditions to 
compute (n) , in contrast to the even lattice case. Indeed, in computing (n) as a function of 
p via eq. (j2I), one finds that (n) the behavior is different for even and odd lattices, ,i.e., (n) 
decays monotonically for odd values of N, while for even values a minimum of (n) for an 
intermediate value p = is observed ^ (see Fig. |21and Tabled] for the detailed analytic 
expressions of (n) as a function of p). In other words, for p > p^^i^ one again observes an 
inverted region where > 0. The value of j9^-„ is strongly size- dependent and shifts to 

the right as N becomes large, i.e., lim^v^ooPmL — ^ (^^e Table |2I). 

A qualitative explanation for the existence of a minimum in the case of an even lattice has 
been given elsewhere js^. This argument suggests that for lattices with an even number of 
sites the number of statistical paths leading to prereactive configurations with both particles 
at nearest neighbor sites is larger than in the odd lattice case. In such a prereactive state, the 
next step will lead to instantaneous reaction with a higher probability if it is asynchronous. 
A small degree of asynchronicity may thus diminish the probability of mutual avoidance 
within the typical interaction zone and thereby result in an increase of the reaction rate 
with respect to the fully synchronous case p = 1. The maximum efficiency will therefore 
correspond to a mixture of synchronous and asynchronous steps. This effect is expected to 
wash out with increasing lattice size, since "diffusional" transport over long distances then 
becomes the rate limiting step, and the latter is optimized by purely synchronous transport. 
Note also that a very similar argument may be invoked to explain the shift of the minimum 
towards larger p-values observed in the curves (n)j(p) with increasing initial separation i. 

These results can be straightforwardly extended to higher-order moments. In the ID case, 
this is most naturally done via a generating function approach [32, .32], but the computation 
of the second-order moment is also possible in the framework of Markov theory (cf. eqs. 
([7j) and 0). Remarkably, the variances {v)i and (v) behave in a similar way as the first 
order moments. For example, the global variance (v) decreases monotonically with p for 
odd lattices, whereas for even values of the curves {v){p) display a minimum at a value 
p = p^^in which again rapidly shifts towards one with increasing lattice size (see Tables El 



This has been verified for lattice sizes up to N=36. 
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and in and Fig. HJ. One has Pmin < Pmin with the exception of the = 4 lattice, for which 
both probabihties are equal (compare Tables |21 and |3]) . Thus, the quantities (f) and (ji) 
cannot be minimized simultaneously in general. 

In the limiting cases p = (where the heavy walker plays the role of a stationary trap) 
and p = 1 (limit of identical walkers), it is possible to recover previous results valid for 
arbitrary lattice size N. In the framework of the Markov approach, these expressions for the 
reaction time and its variance can be constructed by identifying patterns in the (integer) 
values generated in calculating these quantities as a function of increasing A^. In particular, 
one has 

(n), = t{N-i) (10) 



in the p = case [3j| and 

{nh = + f + - (-1)^)(-1^ + \ (iV + ^^^^) (1 - (-in (11) 

in the p = 1 case, whereby i can in both cases take any possible values of the initial 
interparticle distance d^, i.e. i = 1,2, . . . , [A^/2]. 

Note that when and i take even values, one has {n)i = i {N — i)/2, implying that the 
reaction time takes half the value of the p = case, as one would expect from the result 
obtained in the continuum limit. If one now considers the average over a uniform ensemble 
of non-reactive configurations one respectively has 

(n) = J2 ^^H^ = 6 (12) 

i=l 

for p = H and 



(n) = A^(iV + l)(A^ + 2)/(12(iV- 1)) even, 

= (A^ + l)(A^ + 3)/12 AT odd (13) 

for p = 1. According to Eqs. ()12j) and ()13p the purely synchronous case {p = 1) is more 
effective than the purely asynchronous case (p = 0) for all lattices with A^ > 3, as one would 
anticipate from the continuum approximation. 
Turning now to the variance, one has 

, , 2N ^ 2. A^(A^2_2) 



and 



^ N{N + l){N-2){N + 2) 
^ ' 30 ■ ^ ' 



when p = 0. In the p = 1 case the local variance is given by 
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(n^)i = —i^ - + (i - \n + 1)(1 - + (—N^ + —N + -Ar^)^ + - + -N 

^ ' 12 6 ^6 4^ A V v^2 12 4^83 

+AAr3 + -N^ _ (1 + iiv + ^Ar3 + iAr2)(-l)' - -{N + l)Ari(-l)^ (16) 
24 4 ^8 3 24 ^ 4 ^ 4^ ^ ^ ^ ^ ^ > 

for even values of and 



-Ln + — + Iat^ + {—N - — - -N^)(-iY + (Iat - — )2(-l)' (17) 

24 24 8 ^24 24 8 ^ M 12^ ^ ^ ^ ^ 

for odd values of A^. Thus, the global variance becomes 

/ \ = / + + + 2Ar + 2)/(120(Ar - 1)) for N even, , 

\ (iV + l)(A^ + 3)(iV2 + 2iV-5)/120 for iV odd. ^ > 



3.2. Confining boundary conditions 



We now consider the behavior of a ID system subject to confining boundary conditions, 
and compare this with the behavior summarized above for periodic boundary conditions. 
The differences can already be seen by examining the three simplest cases: = 4, = 5 
and A^ = 6. 

In the N = A case there are 6-symmetry distinct initial states in the A^ = 4 case, as 
already seen in 12.11 The initial states 1, 2 and 3 correspond to configurations with walker 
2 at a boundary site. In all three cases, the associated reaction times are found to decay 
monotonically with p. In contrast, the {n)i{p) plots for the initial conditions 4, 5 and 6 
display at least one inverted region where ^^j^ > (see Fig. E}. The behavior for initial 
condition 4 is rather complex, i.e., the reaction time displays a maximum for small values 
of p and a minimum for large values of p, thereby giving rise to a double inverted region, 
as shown in the inset of Fig. On the other hand, (n)^ rapidly reaches a maximum for 
small values of p and then decays sharply with p. Finally, {n)^ increases monotonically for 
all p, i.e., the inverted region comprises the whole p-interval. We thus see that the behavior 
of the efficiency in the confining case is far more complex than in the periodic case, where 
only monotonically decreasing curves or curves with a single minimum are observed. 

In the cases A^ = 5, 6 a similar qualitative behavior is observed. More specifically, the 
initial states 1 to A^ — 1 correspond to monotonically decreasing reaction times, while the 
initial states A^ to A^s are associated with {n)i{p) curves displaying inverted regions with 
positive p-derivative. In some cases the inverted regions extend over a very short p interval, 
this is e.g. the case for the initial states 6 and 7 in the N = 5 lattice and the initial states 
7, 8 and 9 in the A^ = 6 case, whose associated (n)i's display a single maximum very close 
to p = 0. The other reaction times displaying anomalous behavior are plotted in Fig. IHl 
It is noteworthy that those initial conditions for which a double inverted region is observed 
correspond to contiguous walker positions, whereby the heavy walker is placed at an interior 
site closer to the boundary or at most at the same distance from the boundary as the light 
walker. In contrast, those configurations with contiguous particles in which the light walker 
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is closer to the boundary are associated with monotonically increasing {n)i{p) profiles, as 
one might intuitively expect (for large p the heavy walker can move away from the light 
walker more easily leaving a larger space between the boundary and itself available for the 
motion of the light walker). Finally, all the other initial conditions giving rise to anomalous 
behavior correspond to (n)j(p) profiles displaying a single maximum as a function of p. 

Let us now compare the above behavior of the local reaction time with our findings for 
the N = 4 periodic case. To begin with, note that if one switches from confining to periodic 
boundaries. States 2 and 5 become equivalent to, say. State 2 in the periodic case. Fig. [7| 
(left) displays a comparative plot for the {n)i{p) curves in these cases. It is seen that for all 
values of p the reaction times {11)2 and {n)^ are well below the corresponding values of (71)2 
for a periodic lattice. On the other hand, states 1, 3, 4 and 6 become equivalent to state 
1 when switching to periodic boundary conditions. As seen in Fig. [7| (right) the reaction 
times (71)1,3,4 in the confining system are larger than the corresponding reaction time (?7.)i 
in the confining case for most values of p. Exceptions are found for {n)i and (n)4, which 
become smaller in the high-p regime. In contrast, the initial state 6 is the only one in which 
for any "temperature" the reaction is more efficient in the confining than in the periodic 
case. 

In the cases = 5, 6, one has a similar behavior of {n)i and {n)^^, i.e., there is a signif- 
icant p interval where (n)i becomes smaller than in the periodic case, whereas this interval 
extends to the full physical range [0,1] if one starts from the A^^-th initial configuration, in 
which the light walker is placed at a boundary site next to the heavy walker. 

Finally, we turn to the difference in behavior between the confining and the periodic 
case at the level of the global reaction time. In the confining case, the weights Wi used 
to compute (n) take a very simple form, i.e., Wi = = 2/N{N — 1). One thus has 

^'^) ~ N{N~i) ^i^=i~^^^'^{''^)i- -'■^ addition, one can also define other global averages by 
averaging over subsets of initial conditions, e.g. 

1 ^ 

(^y'^'=w^i E (-)(^?.) (19) 

is the average over a homogeneous subset of initial conditions with a fixed value of the heavy 
walker's initial position. (In the periodic case these averages are clearly z-independent and 
equal to the all averaged reaction time (n)). In the case of a A^ = 4 lattice, there are 
only two such symmetry-distinct averages, namely {n)^'^^^ = (n)^^ and (n)^'2"^ = {n)^^ 
with the heavy walker respectively placed at a boundary site (BS) and an interior site (IS) 
initially. For an arbitrary lattice size A^, these global reaction times can be straightforwardly 
calculated from the local reaction times for each symmetry-distinct initial condition by a 
suitable weighting of the individual contributions. In the (particularly simple) case of the 
A^ = 4 lattice, all weights are equal, i.e.. 



{n)^' = + + (20) 
3 

3 

As it turns out (n)^'^ decreases monotonically with p, while (n)^^ first increases and then 
decreases (see Fig. IS}. Notice that the inequality (n)^"^ > {n)^^ holds for all values of 
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p. In addition, one finds that {n)^^,{n)^^ and thus also (n) — ^ — — are larger than 
the global reaction time (n) for the periodic case, thereby emphasizing that at this level of 
description the confining system is less efficient over the whole temperature range. In the 
cases iV = 5, 6 a similar behavior of the global reaction times is found, i.e., for a given p the 
global reaction times decrease when the heavy walker is shifted to the interior of the lattice, 
but in all cases they remain larger than in the periodic case. 

As for the p-dependence of the global reaction times in the = 5, 6 cases, the reaction 
time {n)^^ = is found to decrease monotonically with p, whereas the other 

global reaction times (corresponding to configurations where the heavy walker starts from 
an interior site) first increase with p, thereupon reach a maximum for p well below 0.5 
and finally decrease for larger p values. These maxima are again lost in the "all averaged" 
reaction time (n) , which decreases monotonically with p for all lattices in the confining case, 
in contrast with the even-odd effect reported for periodic lattices. 

4. 2D LATTICE 

4.1. Periodic boundary conditions 

As in the ID case, one can take advantage of symmetry to simplify the formulation of 
the problem, here by lumping the A''(A'' — 1) possible non-absorbing initial states (a;?,^;^) 
in a square planar lattice into a set of A^^ states. Each symmetry-distinct initial state is 
now characterized by the horizontal distance and the vertical distance dy between the 
walkers (measured in lattice spacings) along the two orthogonal directions x and y specified 
by the vector basis of the elementary lattice cell. Both and dy may take values from 

to [\/N/2], but the absorbing state with both walkers at the same site {d^ = dy = 0) 
is excluded. Taking advantage of exchange symmetry with respect to d° and o?°, one can 
therefore also exclude configurations for which, say, dy > c?°. Thus, the total number of 
symmetry-distinct states is ^ A^, = CR^^^^^^^ - 1 = (1^/21+2) 

— 1. To remain consistent 

with the numbering criterion used in the ID case, we shall number these states from 1 to 
As by increasing Euclidean distance d'^ — {d^Y + {d^Y- For example, for a 3 x 3 lattice 
one has 2 symmetry-distinct states, namely: 

State 1 State 2 

2-1-0 2-0-0 
0-0-0 0-1-0 

0-0-0 0-0-0 

In order to compute the reaction time and its variance, one now proceeds as in the ID case, 
i.e., one first computes the corresponding A^ x A^ fundamental matrix and then respectively 



* In the following equation CR^ stands for the number of combinations with repetition of A elements 
choose B. 
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applies the equations © and (|7j) for the local quantities (n)j and {v)i and eqs. ^ and Q 
for the global averages (n) and (v), whereby the weights Wi are again straightforward to 
compute. 

Let us first discuss the behavior of the individual contributions (n)j to the global reaction 
time (n). In the case of a 2 x 2 lattice, there are only two possible symmetry-distinct initial 
states, i.e.. State 1 is a configuration in which the walkers are placed at nearest neighbor 
sites, whereas in State 2 is they are placed diagonal to each other. Fig. El (left) displays the 
behavior of the reaction times {n)i and {71)2 as a function of p. While (n)i decreases with 
increasing p, {11)2 increases with p at a slow rate for low values of p and more rapidly for 
larger values. In the case of a 3 x 3 lattice, the behavior of the reaction time is also different 
for each of the two symmetry-distinct initial conditions 1 and 2, i.e., {n)i displays a shallow 
minimum as a function of p, while {n)2 is a monotonically decreasing function (see Fig. El 
right). 

A careful examination of the behavior for somewhat larger periodic square planar lattices 
(16 < N < 36) suggests the following rule of thumb: those initial states for which the 
effective distance = d}^+dy is even correspond to monotonically decreasing {n)i{p) profiles, 

while the profiles {n)i{p) associated with odd values of (i° show an inverted region. Each of 
these inverted regions is associated with a local minimum with respect to p, except in the 
2x2 lattice case, where the inverted region extends over the whole p interval. As expected, 
the minimum shifts to the right with increasing value of (see Figs. ^land^Jfor the cases 
of a 4 X 4, 5 X 5 and 6x6 torus). This even-odd behavior in terms of d^ is analogous to the 
even-odd behavior in terms of (i° observed in ID (cf. Fig. 

Departures from the above rule become more frequent with increasing lattice size. An 
example is given by the initial state 5 corresponding to the 5x5 square planar lattice, which 
displays a tiny region with ^ > and a very shallow minimum near p = 1, as a careful 
inspection of the thick solid line in the right graph of Fig. ^1 would reveal. Shortly we shall 
see that as a result of this inverted regions at the level of (n) also arise for sufficiently large 
odd values of A^, in contrast with the ID case. 

Next, let us analyze the p-behavior of (n) given in Table El To this end, we first focusing 
on the case of the smallest possible lattice (2x2 lattice, corresponding to A^ = 4). Fig. ^1 
displays the behavior of the global reaction time (n) as a function of p for a 2 x 2 periodic 
square lattice. The inset shows that there is a small region of p- values where < 0. 
However, if one compares this plot to the A^ = 4 plot depicted in Fig. HI (left), it is seen 
that the inverted region with > is much broader in the case of a 2 x 2 lattice than in 

the case of a A^ = 4 ring, leading to a smaller value of p^J^l^ in the 2D case. This behavior 
persists if one considers larger lattices with an even number of sites (compare e.g. the values 
of p^ln for A^ = 16 and A^ = 36 in Tables |2 and IHl) • We thus conclude that for fixed A^ 
the inverted region 1] observed for even lattices in ID is considerably extended in 2D. 

Moreover, in contrast to the ID case an inverted region also appears in the 5x5 lattice 
case (see the corresponding plot in Fig. IT^ . and numerical evidence suggests that inverted 
regions also arise in the 7x7 and the 9x9 cases [23]. 

There is another respect in which the anomalous p-behavior of (n) is here enhanced with 
respect to the ID case. Let us consider the relative enhancement of the efficiency A in the 
case of an even lattice, defined as follows: 
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. ^ 1 [' d{n) {nKp=l)-{n){p = pil) 

{n){p = ptl)Jvtl dp {n){p = p^±) 

From Table [7| it is seen that this quantity is very large with respect to the ID case, even 
for large lattices (A^ > 16). On the other hand, the interval [Pm]„, 1] corresponding to the 
inverted region contracts dramatically with increasing lattice size, as seen from Table IHl and 
expected from the arguments for the 2D case. Thus, a minute amount of asynchronicity 
5p < Q brings about a drastic decrease in (n) with respect to the fully synchronous case 
p = 1. 

We have also studied the global variance (f ) in 2D and find a similar enhancement of 
the inverted regions with respect to the ID case. The corresponding expressions for (w) 
and p^^in listed in Tables |H1 and IHl The values of are very close to those for 
with the exception of the 2x2 lattice, for which one respectively has p^^in ~ 0.1847 and 
pit ^ 0.1042. 

4.2. Confining boundary conditions 

In order to complete the picture of the confining case, we have investigated the reaction 
efficiency for two cases, namely the 3x3 and the 5x5 lattice. 

Let us first discuss the 3x3 lattice. There are 12 symmetry-distinct states altogether, 
i.e., 5 states with the heavy walker at a corner site CS: 

State 1 State 2 State 3 State 4 State 5 

2-1-0 2-0-0 2-0-1 2-0-0 2-0-0 

0-0-0 0-1-0 0-0-0 0-0-1 0-0-0 

0- 0-0 0-0-0 0-0-0 0-0-0 0-0-1 

According to the numbering used in Fig. [T] (right), the above states would correspond to 
those for which X2 = 1. There are 5 more states in which the heavy walker is positioned on 
a midpoint boundary site MBS. These states with X2 = 2 can be pictorially represented as 
follows: 

State 6 State 7 State 8 State 9 State 10 

1- 2-0 0-2-0 0-2-0 0-2-0 0-2-0 
0-0-0 1-0-0 0-1-0 0-0-0 0-0-0 
0-0-0 0-0-0 0-0-0 1-0-0 0-1-0 

Finally, there are 2 more states in which the heavy walker sits in the centrosymmetric site 
CES (i.e., one has X2 = 5), namely 
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State 11 



State 12 



0-0-0 0-0-0 
0-2-0 0-2-0 

0-0-1 0-1-0 

The p behavior of the local reaction times is depicted in Fig. El The reaction times 
associated with the CS initial configurations 1-5 turn out to be monotonically decreasing in 
p. In contrast, the MBS configurations 6-10 and the CES configurations 11-12 are associated 
with reaction times displaying a single maximum. In the MBS case the maximum is located 
at relatively small values of p, whereas in the CES case it is found at somewhat larger p values. 
Interestingly, there are no initial conditions which give rise to a double inverted region, in 
contrast with the ID case. These results are also confirmed by numerical simulations (data 
not shown). 

We now compare the above results with the periodic case. If one connects the boundary 
sites of the confining lattice in a periodic fashion, one finds that States 1, 3, 6, 8, 10 and 12 
collapse into State 1, i.e., the state in which the walkers are positioned aside of each other in 
the periodic lattice. On the other hand. States 2, 4, 5, 7, 9 and 11 collapse into State 2, i.e., 
the state in which both walkers are positioned diagonal to each other in the periodic lattice. 
The plots corresponding to the reaction times in the periodic case are displayed in Fig. HH 
and compared to the corresponding plots for the reaction times in the confining case. For 
all values of p, the local reaction times are smaller in the periodic case than in all confining 
cases ^. This behavior clearly differs from the ID case, where for certain initial conditions 
and p-intervals the efficiency is larger in the confining case. 

Finally, we consider the behavior of the global averages. Consonant with what was done in 
13.21 one can define via equation (fTUl Hhree global averages (n)^^ = {n)^'2^^,{n)'^^^^ = {n)^'^^'^ 
and {n)'~^^ = {n)^'^^^ for all three symmetry-distinct positions of the heavy walker. In the 
present case, these averages are related to the local reaction times as follows: 

cs ^ 2(n)i + {n)2 + 2(^)3 + 2{n)^ + {n)^ 
2(n)6 + 2(n)7 + {n)^ + 2(71)9 + {n)iQ 



^ ' 8 
(^\CES ^ (ri)ii + (n)i2 
^ ' 2 ■ 

The reaction time averaged over all initial conditions is 



(23) 



^ with the exception of two initial configurations for which in the p ~ hmit the reactions times are found 
to be equal to those for the corresponding periodic cases 
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Fig. ^Jdisplays the p-behavior of the three global averages and (n). As in the ID case, the 
efficiency of the reaction for a given p increases as the heavy walker is moved away from the 
centrosymmetric site. On the other hand, one can see that (n)*^"^ decreases monotonically 
with p, whereas (n)*^^*^ and {n)'~^^^ display a maximum. These maxima are averaged out 
when Eq. (j^H) is used to compute (n), leading as in the ID case to monotonically decreasing 
behavior of (n) . The global averages are all strictly larger than in the 2D periodic case, with 
the exception of {n)'"^^ in the limit p — (see Fig. n3|) . 

We have also examined the behavior of the global averages in the 5x5 case. In this 
case there are 6 global averages corresponding to the 6 possible symmetry- distinct initial 
positions of the heavy walker characterized by Euclidian distances d!(j^g = 0, 1, ^/2, 2, v^, 
from the centrosymmetric site. The results obtained are entirely consistent with the findings 
for the 3x3 lattice, i.e., the efficiency increases when the heavy walker is moved towards 
the interior of the lattice. For d^^g < V2 the curves {n)i(p) display a a maximum in p, 
whereas for larger d^^g the decay is monotonic. Once again, the reaction time scaled by 
(n) and computed by averaging over all initial conditions, is monotonically decreasing. 

5. GENERAL CASE 

In this section we shall discuss how results for the periodic lattice case are modified if 
one allows for the possibility that walker 1 does not jump at a given time step. To this end, 
let us now assume that walker 1 and walker 2 respectively jump with probabilities pi and 
P2 at each tick of the clock. Results for this case can be easily obtained using the approach 
developed in 12.21 Note that the initial scenario is recovered as a special case if one sets 
Pi = I and p2 = p. 

5.1. ID case 

Let us ffist examine the behavior of the local reaction times {n)i{pi,p2) as a function of 
P2 when the value of pi is decreased below one. As our previous results for the pi = I case 
show, there is a monotonic decrease in the reaction time with respect to p2 when the state 
number i takes even values (cf. Fig. ^ and minima in p2 for odd values of i. Now, if one 
now gradually decreases pi, the qualitative behavior of the {n)i{p2) proffies does not change 
for even values of i, i.e. one still has < for all values of p2. In contrast, there is a 

qualitative change in behavior for even values of i; the inverted regions with > vanish 
and one finds monotonically-decreasing proffies in this case (see Fig. ^[ left). The even-odd 
effect in i is thus lost. 

Let us now examine the p2 behavior of the global reaction time (n) when pi is decreased. 
Analytic expressions for (n) as a function of pi and p2 for periodic lattices of increasing 
size are given in Appendix A (see Table ITUj) . These expressions are symmetric with respect 
to exchange of pi and p2- Here, for odd values of (n) decreases monotonically with p2 
regardless of the value of pi, i.e., < for all values of p2- In contrast, the behavior in 
the case of an even lattice again depends on the value of pi. As we have seen, for pi = 1 one 
has an inverted region where > 0. However, if one gradually decreases pi the minimum 

of {n){p2) is shifted to the right and below a characteristic A^-dependent threshold value p^fl 
the inverted region finally vanishes and the curves {n){p2) become monotonically decreasing. 
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The critical values p^ l are displayed in Table ^2 for periodic lattices of increasing size. 
These values are obtained as follows. Denote by p'^l^in the value of p2 for which a minimum 
of (n) as a function of p2 is observed. The value of p^2niin a function of pi is then 

(n) 

easily obtained from the condition = for p2 = P^'Jnin- has '^g^™'" < 0,i.e., the 

probability p^2m,in shifted to larger values with decreasing pi until it becomes equal to one. 
This condition yields a polynomial equation in pi with a single root in the interval [0,1]; 
this root corresponds to the threshold value pfl below which the profiles {n)i{p2) decrease 

(n) 

monotonically. Note that one has p\ ^ ^ 1 as N ^ oo. 

Fig. El (right) shows the suppression of the inverted zone in the case of the = 4 lattice. 
Several plots of {n){p2) corresponding to decreasing values of pi are displayed. It is seen that 
the minimum of the curves shifts to the right with decreasing pi until it eventually vanishes 
for pi < p^^l ^ 0.834 (cf. Table ITT|) . Interestingly, the reaction time decreases with pi for 

a fixed value of p2 in the region where p2 ~ 1, while the opposite is true for not too large 
values of p2- This behavior is also confirmed by numerical simulations (data not shown). 

A similar suppresion of the inverted region with decreasing pi is also observed at the level 
of the global variance (v). The corresponding polynomial expressions as well as the critical 
values p^^l are listed in Tables [T21 [131 and jTH As it turns out, the value of p^il is almost 
identical to (see Tables ITT1 and ITU) . 

5.2. 2D case 

As in the ID case, inverted p2 regions are suppressed with decreasing pi, both at the level 
of the local and the global reaction times. Fig. El displays the behavior of the (n)i and (n) 
as a function of p2 for the 2x2 lattice. In both cases, the curves become monotonically 
decreasing below a given threshold value of pi. In contrast, the curves {n)2{p2) remain 
monotonically decreasing for any values of pi. The p-behavior of (n) is thus qualitatively 
similar to the ID case (cf Figs. ITHI and IT7|) . The analytic expressions for (n) in terms of 
Pi and p2 as well as values of p['^l for even square lattices of increasing size are listed in 
Appendix B (Tables UHl and Table ITH|) . 

For 2D lattices the limiting behavior — > 1 as — > oo is the same as for ID. Here, 

however, for a given number of lattice sites the values of displayed in Table El are 
systematically smaller than those shown in Table ^2 This fact suggests that in higher 
dimensions the inversion effect is more robust vis a vis a decrease in the mobility of walker 
1. Finally, suppression of the inverted region with decreasing pi is also seen in the behavior 
of the global variance (v). The polynomial expressions for (v) and the associated values of 
p^il are displayed in Tables IT7I and ITHl 

6. DISCUSSION 

Let us now summarize the main results obtained so far and introduce some physical 
arguments to explain them. 
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6.1. Results for periodic lattices 

The existence of an inverted region in the case of a ID periodic lattice is governed by 
a parity effect. For even values of the interparticle distance dP no inverted regions are 
observed. However, for odd values of dP the reaction time as a function of p displays a single 
minimum for a p-value larger than 1/2, which rapidly shifts towards one with increasing 
and increasing lattice size. In 2D, a similar even-odd effect holds only for sufficiently small 
lattices; for larger lattices exceptions to this rule become increasingly frequent. On the other 
hand, the inverted regions become broader with respect to the ID case and the enhancement 
of the reaction's efficiency within the inverted region also becomes much larger (cf. Table 

El). 

As already anticipated in 13.11 the above even-odd effect in ID is closely related to the 
fact that in prereactive configurations where both walkers are located at contiguous sites an 
asynchronous step leads to reaction with a probability twice as large as for a synchronous 
step (1/2 vs. 1/4). For such a configuration asynchronous dynamics due to a lower system 
temperature minimizes the probability of mutual avoidance once the particles are within the 
interaction zone (in which the next time step may lead to reaction). For sufficiently high 
values of p the number of statistical paths leading to such nearest-neighbor configurations 
are higher for odd values of dP than for even values. In this high p regime, the reaction time 
increases with p for odd , while it decreases with p for even values of d^. The number of 
initial conditions with odd values of d^ is larger for an even value than for an odd value of 
N] thus the total increase in p of the associated reaction times in the former case overcomes 
the p-decrease of the reaction times associated with even-valued initial distances and the 
global reaction time (n) increases. 

In ID, the size of the inverted region for odd values of the interparticle distance d^ must 
decrease with increasing d^ (Fig. 12)) or increasing lattice size N (Fig. E)), since discrete 
diffusion over long distances becomes the rate limiting step. In 2D, the enhancement of the 
even-odd effect with respect to the ID case can be understood as follows. For a given number 
of sites iV, the initial-condition-averaged distance between the particles is smaller than in ID, 
thereby reducing the role of long-distance hopping with respect to mutual avoidance in the 
interaction zone. In addition, the probability of reaction for nearest-neighbor configurations 
in a sufficiently large lattice is four times higher for an asynchronous than for a synchronous 
step (1/4 vs. 1/16). 

Finally, we have seen that the inverted regions in the periodic case are suppressed if the 
mobility of the light walker lies below a size-dependent threshold value. This is reasonable, 
since a prerequisite for the onset of the even-odd effect is that both particles must first have a 
minimal mobility to come into the vicinity of each other. If this is not the case, the dynamics 
is controlled by the relative diffusion coefficient, which is a monotonically increasing function 
of pi and p2. 

6.2. Results for confining lattices 

In small confining lattices the behavior of the reaction time is even more complex. In 
ID, the reaction time may decrease monotonically, increase monotonically, display a single 
maximum or even a maximum and a minimum, giving rise to inverted regions in the last 
three cases (cf. Fig. Those initial conditions with the heavy walker at a boundary site do 
not give rise to inverted regions, whereas all the others do. Initial configurations where the 
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light walker is on a boundary site and the heavy walker is next to it result in reaction times 
which increase monotonically with p, thereby giving rise to an inverted region of maximum 
size. In particular, the p = case (one walker plus trap case) then becomes more efficient 
than the p = 1 case (case of identical walkers). Double inverted regions are observed for 
those initial states where both walkers are at interior sites next to each other with the light 
walker closer to the lattice boundary or at the same distance from the boundary as the 
heavy walker. In contrast, no double inverted regions are observed in the 2D cases studied. 

Both in the ID and the 2D case, the global reaction time (n)^2=* (average over the 
positions of the light walker) displays a local maximum as a function of p as long as the 
heavy walker starts sufficiently close to the geometric center of the lattice, and it decays 
monotonically with p as the heavy walker is shifted away from the center (cf. Figs. |H1 and 
[T3j) . The reaction time (n) is found to decrease monotonically with increasing p, in contrast 
with the case of a periodic lattice with an even number of sites. 

This complex behavior for confining boundary conditions can, at least in ID, be under- 
stood in terms of an description based on the concept of a "fluctuating" lattice. To illustrate 
the use of such a metaphor, consider the ID case subject to the choice p = 0. Here, walker 
2 plays the role of an immobile trap which divides the lattice into two disconnected sublat- 
tices. Depending on its initial position, walker 1 will perform a random walk between the 
reflecting and the absorbing boundary of one of these sublattices. From the standpoint of 
the reaction time (see Fig. ITHj) . the motion of walker 1 in each sublattice is equivalent to 
its motion in a periodic lattice. For most initial conditions, the size of the corresponding 
periodic lattice is larger than the size of the original confining lattice, implying that for a 
given value of the reaction time (a monotonically increasing function of N) also becomes 
larger than in the periodic case. Thus, in the small p limit the reaction is less efficient in 
the confining case than in the periodic case. An exception is found for the N{N — l)/2-th 
initial condition (see e.g. the curve for {n)^ in Fig. [7[ right), for which the effective size of 
the associated periodic lattice is smaller than A^. Such configurations are no longer found 
in the 2D cases we have examined, where periodic lattices are more efficient than confining 
ones, with the exception of two cases in the limit p ^ for which the efficiency becomes 
the same (cf. Fig. HD). This behavior is due to the fact that the trapping effect exerted 
by the heavy walker and the confining lattice boundary on the light walker is lost in higher 
dimensions. 

If one now slowly increases p, the above effective description in terms of a periodic lattice 
remains valid, since walker 2 remains a deep trap. However, the fact that walker 2 can 
now move leads to time fluctuations of the lattice size of the associated periodic lattice. In 
the case where walker 2 starts off from a boundary site (corresponding to the ffist A^ initial 
conditions), fluctuations can only lead to a smaller effective size than the initial one, resulting 
in a monotonic decrease of the reaction time with respect to the p = case. However, for 
any other initial condition with walker 2 at an interior site small fluctuations in the position 
of the trap can in many cases lead to a larger effective size, resulting in an increase of the 
reaction time. On the other hand, for somewhat larger values of p excursions of walker 
2 towards the boundaries become more frequent, resulting in a more efficient confinement 
of walker 1 and thereby in a decrease of the reaction time. In other words, the effect due 
to lattice size fluctuations is overcome by the significant increase in the relative diffusion 
coefficient, which is a monotonically increasing function of p [21]. The reaction time, which 
in a continuum description is inversely proportional to the diffusion coefficient, will therefore 
decrease with p in this regime. As a result of the combination of the above low-p and high-p 
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effects, for those configurations where walker 2 is sufficiently far from the boundary and the 
initial particle separation is sufficiently large (such that the continuum picture applies for 
sufficiently large p) the reaction time first increases and then decreases with p. 

The additional increase of the reaction time with p in the high p region for the subset of 
initial conditions associated with a double inverted region can possibly be explained in terms 
of an enhancement of mutual avoidance effects when particles are close to each other and 
far from the boundary ^, resulting in an increased escape probability. The absence of this 
additional high-p inverted region in the 2D case is probably due to the existence of additional 
degrees of freedom which no longer make a description in terms of the periodic case possible. 
One could therefore expect that the even-odd effect leading to the high-p inverted region 
does no longer hold and that the dynamics is essentially governed by long-range diffusion, 
resulting in a p-decrease of the reaction time. 

In ID, the disappearance of the inverted region associated with the the global encounter 
time (?T,)^2=« as walker 2 is shifted away from the inmost lattice sites to the boundaries could 
be explained by the fact that lattice size fluctuations responsible for the low p inverted 
region play a less important role when the initial distance between the walkers becomes 
larger. Clearly, the average of the interparticle distance over the positions of walker 1 
becomes larger when walker 2 is shifted away from the centrosymmetric site. 

Finally, the monotonic decrease of the global reaction time (n) over the whole p range can 
be explained as follows. In the small p limit the effective shortening of the lattice through the 
motion of walker 2 from the boundary to the bulk for the first N initial conditions overcomes 
the effective increase in system size for the other initial conditions. In the opposite limit 
of large p, the even-odd effect leading to mutual avoidance in the effective periodic-lattice 
description is weakened due to the lattice size fluctuations, and the global reaction time 
decreases. 



7. CONCLUSIONS 

The principal conclusion that can be drawn from the study presented in this paper is that, 
in characterizing the efficiency of (two-channel) reactions between independently-mobile 
reactants on a lattice, for small systems one may find regimes where, counter one's intuition, 
the reaction efficiency can decrease with increasing "temperature." Moreover, both minima 
and maxima in the reaction efficiency have been found and quantified, thus establishing 
the existence of one or more "inverted regions." We stress that, in contrast to Marcus 
theory, the emergence of these anomalous regions is a consequence of a (strictly) classical 
interplay between system geometry and individual particle dynamics (and the importance 
of interactions within a certain reaction radius or "zone"). We have demonstrated and 
quantified that these inverted regions arise in ID and 2D lattices for a wide class of initial 
conditions. Also, we note that numerical evidence has been presented in |23| that the effect 
in present in 3D, and for cases where the reaction "rules" are different from the ones adopted 
in this work 

So that the plethora of results reported here not obscure a principal insight, we emphasize 
that the discrete character of the geometric support is an important (perhaps critical) ingre- 
dient in generating the emergence of an inverted region. It is precisely the discrete geometry 



^ Recall that double inverted regions are associated with initial conditions which fulfil these two conditions. 
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of the lattice which allows one to distinguish between reaction channels (viz., irreversible 
termination of the reaction owing to same site occupancy or nearest-neighbor crossing). 
Such distinctions are lost in the continuum limit, where the efficiency of a diffusion- reaction 
process is characterized universally by an increase in reaction efficiency with respect to an 
increase in system temperature. 

In the cases we have studied here, much of the observed anomalous behavior can be un- 
derstood in terms of a competition between diffusional transport over long distances and the 
probability of mutual avoidance inside the interaction zone. A very large jump probability 
p may favor transport over long distances, since it increases the relative diffusion coefficient. 
However, it may also decrease the residence time in the interaction zone. As we have seen, 
the latter effect may become dominant in sufficiently small lattices and thus determine the 
anomalous decrease of the efficiency with temperature. Remarkably, the fluctuation-driven 
resonance and antiresonance phenomena reported in this work do not involve coloured noise, 
in contrast with other first-passage problems such as resonant activation [s^. 

The study presented here was motivated to a large extent by the remarkable and un- 
expected results found in studying experimentally systems at the nanolength scale. It is 
premature to suggest that the effects uncovered here might be found in a particular experi- 
ment, but it is certainly not premature to note that the sophistication of these experiments 
has now reached a point where the presence of an "inversion regime" might be explored. 

From a theoretical point of view, there are many facets of the present problem that need 
to be developed. For example, one can study the effect of modifying the reaction rules. In 
this context, the case of a reaction probability pfi < 1 per collision is particularly interest- 
ing. An enhancement of the anomalous behavior can be expected in this case, since the 
efficiency of recoUisions in the interaction zone can be strongly diminished by amplifying the 
simultaneous motion of the two reactants. Such counterintuitive effects are rather common 
in binary reaction systems. For example, Saxton has recently considered a binary reaction 
with pji < 1 between two simultaneously moving walkers in a percolation cluster with block- 
ing sites js^. This model was recently proposed to explain protein motion in lipid rafts. 
Saxton came to the interesting conclusion that at sufficiently low reaction rates increasing 
the obstacle concentration might result in an anomalous decrease of the reaction time due 
to an enhancement of recoUisions in the interaction zone. Similar effects are currently being 
investigated for binary reactions in other types of disordered media, e.g. so-called small- 
world networks, in which a very small number of long-range connections between sites may 
have drastic effects upon the reaction time. 

Finally, since non-nearest neighbor jumps [23] are likely to become of increasing impor- 
tance with an increase in the system temperature, it is of some interest to incorporate this 
facet into the present model. Moving from a scenario in which one analyzes the problem in 
terms of a symmetric random walk, incorporating " drift" and accounting for deactivation by 
regular or periodic impurities in the system opens up connections to an even broader range 
of experimental phenomena. It may be hoped that the new insights arising from the study 
of the above systems along with non-classical approaches will help to improve the current 
understanding of encounter-controlled processes in small-size systems. 
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Tables in the main body 



N_ (n) 

2 2/(2 -p) 

3 2 

4 (10/3)(3p-4)/(p2 + 2p-4) 

5 4(2p-5)/(p2 -4) 

6 (28/5) - lOp + 10) /(p3 - 4p2 -Ap + 8) 

7 (4/3)(p2 + 8p-14)/(p2-2) 

8 (12/7)(13p3 + - l26p + n2)/{{p - 2){p^ + Qp^ - 8)) 

9 10(2p3 - - 16p + 24)/((p2 + 2p - 4) [p^ - 2j3 - 4)) 

10 (22/9)(7p'' - 76p3 + 16p2 + 288p - 240)/(p5 - - 12p3 + 32p2 + I6p - 32) 



TABLE 1: Analytic expressions for (n) (ID periodic lattices of increasing size N) 
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N 


(n) 
Pmin 


2 


.000000 


4 


.666667 


6 


.859648 


8 


.920393 


10 


.948308 


12 


.963611 


14 


.972951 


16 


.979086 


18 


.983337 


20 


.986407 


22 


.988698 


24 


.99045 


26 


.99183 


28 


.99292 


30 


.99381 


32 


.99455 


34 


.99515 


36 


.99567 



TABLE 2: Values of (ID periodic lattices of increasing size N) 
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N (v) 

"2 2p/{2-pf 

3 2 

4 (2/3) (192 - 316p + 152p2 - 15p3)/(p2 + 2p - 4)^ 

5 4(-2p3 + 25p2 - 80p + 84)/(p2 _ 4)^ 

6 (28/5)(-1096p + 768p2 - 206^^ + 28p^ - / + 512)/(8 - 4p + p3 _ 4p2)2 

7 (4/3)(-320j9 + 88^92 + 8^93 +/ + 252)/(p2 - 2)^ 

8 (12/7)(-13p'^ + 252/ + 546p5 - 4736^1 - 2184^3 

+35520p2 - 50848p + 21504)/((p3 + 6p2 - %f{p - 2f) 

9 2(-10p'^ + 253/ - 816/ - 2748^"^ + 11872/ + 5552/ 

-42496p + 29568)/((p2 -2p- 4)^{p'^ + 2p- 4)^) 



TABLE 3: Analytic expressions for {v) (ID periodic lattices of increasing size N) 



N p^^'l 

fmin 

2 .000000 

4 .666667 

6 .852189 

8 .914315 

10 .943742 

TABLE 4: Values of (ID periodic lattices of increasing size N) 
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N_ (n) 

1 (2/3)(40 - 31p)/(-p2 - 6p + 8) 

9 8(9-2p)/(8-p2) 

16 (8/45)(1191p3 + 252p2 - 14252;? + 13184)/(3/ + 24^^ _ 5Qp^ _ 9Qp + 128) 

25 (76/3)(-9p3 + 2p^ + 320 - 124^^ - 48p)/{p^ + Ap'^ - A4p^ - UAp'^ + 6Ap + 256) 

36 (2/1225)(169521024p3 - 99291648^2 - 21958740/ + 760284/ 

+247652352 - 296662016p + 494505/ + 614144p'^)/(5376/ 
-8704/ - 1176/ - 256/ - 6144p + 8192 + 54/ + 5/ + 2688/) 



TABLE 5: Polynomial expressions for {n) (2D periodic square planar lattices of increasing size 
N) 



N 




Pmin 


4 




.104208 


16 




.750207 


36 




.840892 


TABLE 6: 


Values of (2D periodic square planar lattice of increasing 


: (even-valued) size A'') 


N 


A (ID) 


A {2D) 


4 


0.111 


0.802 


16 


0.067 


0.416 


36 


0.016 


0.344 



TABLE 7: Relative enhancement A of the global reaction time for ID periodic lattices and 2D 
periodic square planar lattices with the same number of sites N. 
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N_ (v) 

1 -(2/3)(1240p + 31p3 - 560p2 - 768)/(p2 + Qp- 8)^ 

9 -8(288p + 2p3 - 4 V - 584) /(p^ - 8)^ 

16 -(8/135) (195677696J9 + 15943296^'' - 5833248p''^ - 115741696^^ - 511812/ 

-705816/ + 10719/ - 88932352)/(3/ + 24/ - 56/ - 96p + 128)^ 

25 -(4/3)(16064512p - 8024128p^ - 3642368/ + 35447808/ - 2255/ - 48480256 

-813472p5 + 187628/ + 10448/ + 38p'^)/(/ + 4/ - 44/ - 144p2 + 64p + 256)^ 

36 - (2/8575) (1690178540273664p - 9697076228pi-^ + 226056069312pi2 + I7307675p^' 



+1052496354082816p^ - 1589772806520832p''^ - 494222977269760^^ + 435557221152p 
+87781881765888/ + 477281389805568/ - 477098773250048/ - 39787096070 144p' 
-7123662261248^1° - 2654087139328/ - 697746159828992 - 2669903720^^^) 
/(-6144p + 2688/ + 5376^^ - 8704/ + 5/ + 8192 - 1176/ - 256p^ + 54/)^ 



TABLE 8: Polynomial expressions for (v) (2D periodic square planar lattice of increasing size N) 



N 


{v} 
Pmin 


4 


.184684 


16 


.748812 


36 


.840409 



TABLE 9: Values of (2D periodic square planar lattice of increasing (even-valued) size N) 
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Figures 
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FIG. 1: Integer coordinates associated with each site for a ID lattice with N = 6 (left) and a 2D 
lattice with N = 9 (right). 
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FIG. 2: p-behavior of the local reaction times in a N = 9 case (left) and the N = 10 case (right). 
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FIG. 5: p-Behavior of the local reaction times for those initial conditions characterized by the 
onset of an inverted region (ID lattice with N = A) The inset clearly shows the double inverted 
region displayed by {n)^. 
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FIG. 6: p-Behavior of the local reaction times for some of the initial conditions characterized by the 
onset of an inverted region. The left (right) graph corresponds to a ID lattice with N = 5{N = 6). 
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FIG. 7: p-behavior of the local reaction times in the confining and the periodic N = 4 case. The 
initial states 2 and 5 in the confining lattice collapse to the initial state 2 ( corresponding to walkers 

at next to nearest sites in the periodic lattice) when the lattice boundary sites are connected to each 
other (see left graph). The initial states 1,3, 4 o.i^'d 6 in the confining lattice collapse to the initial 
state 1 (corresponding to contiguous walkers in aperiodic lattice) when the lattice boundary sites 
are connected to each other (see right graph). 
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FIG. 8: p-Behavior of the global reaction times for the periodic and the confining N = 4 case. 
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FIG. 9: p-Behavior of the local reaction times (n)i and (71)2 in a 2 x 2 (left) and a 3 x 3 periodic 
square planar lattice (right). The values of the effective distance dP associated with the initial state 
1 and 2 are respectively 1 and 2. 




FIG. 10: p-Behavior of the local reaction times for each of the five symmetry- distinct initial 
configurations in a 4 x 4 square planar lattice (left) and a 5 x 5 square planar lattice (right). In 
both cases the values of d^ associated with the initial states 1, 2, 3, 4 o,nd 5 are respectively 1, 2, 
2, 3 and 4. 
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FIG. 11: p-Behavior of the local reaction times for each the 9 symmetry- distinct initial configura- 
tions m a 6 X 6 square planar lattice. The values of dP associated with the initial states 1, 2, 3, 4, 
5, 6, 7, 8 and 9 are respectively 1, 2, 2, 3, 4, 3, 4, 5 and 6. 




FIG. 12: p-behavior of (n) in an N = A periodic square planar lattice. The inset shows the 
minimum of the curve for p 0.104. 



34 




0.8 1 f 

0.0 0.2 0.4 0.6 0.8 1.0 

P 



FIG. 13: p-Behavior of the global reaction time (divided by the lattice size) for 2D square planar 
lattices of increasing size. 




FIG. 14: p-behavior of the local reaction times in the confining 3x3 case. The initial states 1, 
3, 6, 8, 10 and 12 in the confining case collapse to the same state when the boundary sites are 
connected periodically, i.e., State 1 with contiguous walkers in the periodic lattice (see left graph). 
The initial states 2, 4, 5, 7, 9 and 11 also collapse to the same state, i.e.. State 2 corresponding 
to walkers placed diagonally in the periodic lattice (see right graph). The thick solid lines show the 
p-behavior in the periodic case. 
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FIG. 15: p-Behavior of the global reaction times in the confining 3x3 case. The thick solid line 
corresponds to the behavior of (n) in the periodic 3x3 case. 




FIG. 16: Behavior of the curves (n)i(p2) (left) and (n)(p2) (right) for different values of pi in a 
periodic ID lattice with N = A. 
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FIG. 18: Correspondence between a confining four-site ID system with a deep trap T and periodic 
systems of different sizes for the two possible symmetry- distinct trap configurations. Sites with 
same numbers are symmetry- equivalent. In our model, the heavy walker 2 plays the role of the 
deep trap T in the figure. For p = the system in the left subfigure will be equivalent to a three- 
(five-) site system if walker 1 starts to the left (right) of walker 2. Similarly, the system in the 
right subfigure will be equivalent to a seven-site periodic system. From this picture one (-am, easily 
infer that the motion of the trap in the case p > is associated with lattice size fluctuations in the 
periodic lattice representation. 
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APPENDIX A: TABLES WITH ANALYTIC RESULTS FOR THE GENERAL 
CASE IN ID 



N (n) 

2 2/(2pi + 2p2 - 3piP2) 

3 2/{pi+p2-piP2) 

4 (10/3) (4pi + 4p2 - 7pip2)/{^pj + 8pip2 - Wpjp2 + ^pI - Wpipl + ^pIpI) 

5 (20pi + 20p2 - 28pip2)/{{piP2 - 2pi - 2p2){3piP2 - 2pi - 2p2)) 

6 (28/5)(-20pip2 - lOp? + 30p?p2 - lOpi + SOpip^ - 2lplpl) / {-2Aplp2 
-2Apipl + 56pjpl - 8pl - 8pl + 28p\p2 - 28p\pl + 28pwl - '^^pIpI + 7p?pi) 

7 (-4/3)(21p2p2 _ 36p2p2 + 14pf - Mpipl + 28piP2 + Up?.)/ 
iiPiP2 - 4p?P2 + 2pl - Apipl + 4pip2 + 2^2) (-Pi - P2 + P1P2)) 

8 (12/7)(-231pfpi + 594pfp| - 462pfp2 + 112pf + mAp\pl - 924pfpi 
+336pfp2 - 462pip3 + 336pip^ + ll2pl) / {{Zp\pl - I8p\pl + 24pfp2 - 8pf 

-18pIpI + 4&pIpI - 2ApIp2 + 24pip3 _ 24pipi - 8pl){-2p^ - 2p2 + ?,pw2)) 

9 -10(33j3fpi - 99pfp^ + 88pfp2 - 24pf - 99pfp^ + 176pfp^ - 72pfp2 

+88pip\ - 72pipl - 2Apl)/ {{pjpl - QpIp2 + 4pf - Qpipl + 8pip2 
+4p2)(4p2 + 8pip2 - 10pfp2 + 4pi - lOpipi + ^pIpD) 

10 (-22/9)(429pfpf - 1716pfp^ + 240j9| + 2288pfp| - 1248pip| - 1716pfpi - 1248p|j92 



+240pf + 2288p|p^ - 3744p2p3 + geOpip^ _ 3744^3^2 ^ 950^3^2 + UAOpjpl + 4576pfpf)/ 
(llpfpl - llOpfpl + 308pfpf - 352pIpI + 176pfp2 - 32pf - IWpfp^ + 616pfp^ 
-W56pfpl + 704p|p| - 160pfp2 + 308pfp| - 1056pfp| + 1056pfp| - 320pfp| 
-352j9|j9f + 704pfp| - 320pfpi + 176pip^ - IGOpipj - 32p^) 

TABLE 10: Polynomial expressions for (n) in the general case (ID periodic lattices of increasing 
size N). 
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N 

T" 

6 

8 

10 
12 



(n) 

.833972 
.902207 

.936378 
.955625 
.967419 



TABLE 11: Critical values (ID periodic lattices of increasing size A^). 



AT 



{v) 



2(2 - 2pi - 2p2 + 3piP2)/(-2pi - 2p2 + ^ViP2y 



3 
4 



2(2 - - P2 + PlP2)/{-Pl -P2+ PlV2f 



(2/3)(272pf + 544pip2 - 1200pip^ + 1530pfp^ + 340pfp2 - 450pfp^ 
+340p|j3i - 450j9fpi + 175pfp| - 1200pfp2 - 80pf - 80p^ 
+272p|)/(4pf + Spip2 - 10pfp2 + 4p| - IQpwl + 5pfpi)2 

4(21p3p3 _ np\pl + 68pIp2 - 20pl - 7\pIpI + MSplpl - Zmp\p2 + IMpl + 68p^pi 
-356pipi + 208pip2 - 20pi + IQApl) / {{p^p2 - 2pi - 2p2)^(-2pi - 2p2 + 3pip2)^) 



(28/5)(-3984pip^ + 10056pfp| 
+1498pfj3| - 11592pfp4 _^ 5754p4^4 



80pf + 592p| 



, ____ 11592p|p3 „ i288pfpi 
+1498p^^j3^ - 11592pfp^ + blMplpl - 798pfp| + 520pfp2 - 3984pfp2 + 10056j3fp| 
+592pf + 3552pfp^ + 19072pfp| - 11552p?p| + 2368pipi - 11552pfj9| + 2368pfp2 
+1498pfp| - 798p|p^ + UIpIpI - 80p| + 520p|j3i - 1288p|pf)/(-24pfp2 - 24pip| 
+56p2p2 + 7p3p3 _ 28p3p2 + 28p3pi - 28p2p3 + 28pfp2 - 8p? - 8pi)2 

(4/3)(-1596pip^ + 3384pfpf - 28pf + 280p| - 'i24Qp\pl - 328pfpi + 320pfpi 
-3240pfpt + l2,Up\p\ - Ulp\p\ + 156pfp2 - 1596pfp2 + 3384p|p^ + 280p| 
+1680pfj9i + 6456pfpi - 4648j3fp| + 1120p|pi - 4648pfpi + 1120pfp2 + 320j3fpi 
-lAlp\pl + 2lp\pl - 2Spl + IbQplpi - 2,2SpIpI) / {{pIpI - 4j3fj52 + 2pl 
-Apipl + 4pip2 + 2pi)^(-pi - P2 + PlP2)^) 



TABLE 12: Polynomial expressions for {v) in the general case (ID periodic lattices of size N 
to N = 7). 
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N {v) 

"S (12/7)(349440pfi?| - 1998080pfpi - 1005312pfp| + 2845920pfpi - 1998080^??^^ 

+4230240pfpf - 3982752pfp^ + 139776pfp2 + 349440pfp| + 465920pfp^ - 54864p|p2 

-5A86iplpl + mmp\pl + lososop^p^ - i92mplp% - immp\pl + 2079p^p^ 

-109512pjp| + 64746pjp^ + lOSOSOpJp^ - 109512pfp^ + 236610^?^^ - 867528pfpf 
+746256p|pf - 1419024p^pf + 1505928p|pf - 867528p?p| - 1419024p5p| + 1505928p5pt 
+746256j9fp| + 23296p^ + 15456pjp2 + 15456p^pi - 1792j3j; + 23296pf - 206080j?f j92 
-206080p^Pi + 2845920pfpi - 3982752p|p| + 2805696pfp| + 139776p|pi 
-1005312p5p2 _ I792p7)/((3p3p3 _ 18^3^2 ^ 24pfp2 - 8pf - ISplpl + ASpjpl 

-24pfp2 + 24p3pi - 24pipi - 8pi)2(-2pi - 2p2 + 3pip2)^) 

9 2(3000448pfpf - 1203840pfp^ + 188928pfp2 - 3669728pfp| + 3000448pfp^ - 1203840pfp^ 

+188928pip^ - 3669728pfp| + 2238576pfp| + 629760pf + 472320p|p^ - 2394240pfpi 
+472320j3fp| - 2394240pfp| + 4463872j3f + 825plpl + 783232j9fp| - 1300176pfpi 
+1195368pfp^ - 591060pfp| - 70340pfp| + 76080pfp^ - 46160pfp^ + 35860pfp^ 
-9075pfp^ - 70340pjpf + 76080p^pi - 46160pjp| + 35860pjpi - 591060pfp| + 1195368pfp| 
-1300176pfpf + 136620p?pf - 9075pjpf + 14720pjp2 + 783232pfp^ + 14720pip^ + 31488p| 
+31488pf - 1920pl - 1920p^ - 246144pfp2 - 246144pip|)/((pfp^ - 6pipl + 4p| - 6pfp2 
+8;ji;j2 + 4pf )2(4p| + 8pip2 - WpiP^ + 4p? - Wplp2 + ^plplf) 

10 (22/9)(8687616pfpi - 57738240pfp| + 8687616pfp^ - 57738240p^[p^ + 150474240pfpi 

-97613824pj'p| + 10859520pfp| + 1236378pfp| + 37019488pjp^ + 20033904pfp| 
-7090512p^p| + 155136p| - 20579328pfpi + 32759936p?p^ - 32651520pfpi - 1691136pfp2 
+20033904^5 jj^ - 7090512pfp5 + 7893504pfp| - 20579328pfp| + 32759936pfpi - 32651520pfp| 
-1691136pip| + 7893504pfp| + 4343808pfpi - 34707456pfp^ + 113143296pfp| - 192586240pfp| 
+155136pf - 97613824p|p^ + 45717504pfp| - 11630592pfp^ + 122302400pfp^ - 90295040]??^^ 

+45717504pjpf - 11630592p[p^ + 122302400^^^1 - 192586240pfp^ + 113143296p|p| 
-34707456pf j?^ + 181026560j3fp^ - 90295040j3{'p| + 1241088j5[j52 + 1241088j3ij9^ + 4343808pfp^ 
-7680p? - 7680pf - 1531904p5'p| + 1524160p?p| - 944944p?p| + 346060p?p^ + 970816p?pi 
-mUAplpl + 82176p?p2 - 66066p?pf - 1531904pfp| + 1524160pfpl - 944944pfp| 
+346060pjp| + 970816pfpl - SlTSUpjp^ + 82176pipf - 66066pfp| + 4719p?p^)/(llpfp^ 
-llOpfpl + SOSpIpI - 352pIpI + UGpipl - 32p^ - llOpfp^ + 616p|p| - lObGpfp^ 
+704pIp^ - 160pipj + 308pfp^ - W56pfpl + 1056pfpf - 320pjpl - 352pfp| 
+704pfpl - 320pIpI + 176pIp2 - lG0pjp2 - 32pf )2 



TABLE 13: Polynomial expressions for {v) in the general case ( ID periodic lattices of size N = 7 
to N = 10). 
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N 




4 


.817442 


6 


.891104 


8 


.929001 


10 


.950475 



TABLE 14: Critical values p^], (ID periodic lattices of increasing size). 
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APPENDIX B: TABLES WITH ANALYTIC RESULTS FOR THE GENERAL 
CASE IN 2D 



(n) 

1 -(2/3)(71pip2 - 40p2 - 40pi)/(13pfpi - 22pipl - 22pjp2 + Spj + 16pip2 + Spj) 

9 -8(llpiP2 - 9pi - 9p2)/{7pIpI - 1QpIp2 - IQpwl + Spj + 16pip2 + 8pl) 

16 -(8/45)(26497pfp^ - 68308pfp^ + 53804pfp2 - 13184pf - GSSOSpjpl + 107608pfpi 

-39552j9fp2 + 5380ipipl - 39552pipl - ISlSipj) / {147 pfpj - 66ipfpl + lOOOpfpl 
-608pfp2 + 128pf - 664p?p^ + 2000pfpi - 1824p?p| + 512pfp2 + lOOOpfpf 
-1824pfpi + 768pIpI - GOSpipj + 512pipl + 128pf ) 

25 -(76/3)(1920pfp^ - 3984pfp^ + 1280pfp2 - 1185pfp| - llSSpfpf 

+ 1940p|p^ + 1940pfpf + USOpipl + 3880pfpf - 3984pfp^ - 1328p|p2 
+320p| + 255pfp^ - 1328p|j3i + 320p|)/(95pfp| ~ 684pfp| - 2560pfpi 

+1788j9f pi - 2160pIpI + 1216pfp2 - 256j?f - 684pfp| + 3b7Qp\p\ 
-Ump\pl + 4&Up\pl - 1280pfp2 + 1788pfp| - 6480pfp^ + 7296pfp| 
-2\mplpl + 4864pfp^ - 2bmplpl + 1216pip| - 1280pfpi - 256pi) 

36 -(2/1225)(-12930686208j9lp| - 12930686208pfp^ + 15328776360pfp| 

-38792058624pfp6 _ 5200699392pfp2 _ 5200699392p|p2 ^ 747i0699776pfp| 

+30453427200pfpi - 68813798400pfp| + 30453427200pfp| + 49807133184p5pf 
-38792058624p5p| - 34406899200p?pi + 49807133184pfp| - 34406899200pfp| 
+12181370880pfp| - 68813798400pfpt - 1733566464p5p2 - 1733566464pip| 
+12181370880pfp| + 40604569600pfpf - 8667832320p?pf - 8667832320pfpi 

+298340359p^p^ - 2385516420pjj3^ - 2385516420^?^^ + 7664388180^^^1 + 7664388180pf 
-247652352p^ + 12451783296pfp^ - 6881379840pfp^ + 2030228480pip^ - 247652352p[ 
+12451783296plpi - 6881379840plp^ + 2030228480plj92)/(193480p|pf - 71680pip| 
-530176pfp| + 263680pfp| - 457880pfp| + 633728pfp| - 71680j3fp2 - 2650880pjp| 
-457880j9?j9i - 530176j9^j9i + 263680j9^j9i + 633728j9?j9| - 2650880pfp^ + 3802368pfp| 
-5301760p?p^ - 5301760pfpf + 5273600pfp| - 2508800pfp^ + 458752pfpi 
-2508800j9f pi + 458752p?pi + 3955200pfp| - 1505280p?pf + 3955200j9fpi 
-1505280pfp| + 229376p?pi + 229376p?p^ + 573440pfpf + 193480pfpi + 8192p^ 
+386960pIp^ - 1373640j9lpi - 1373640p5p^ + 2534912pjpi + 2534912pfp^ 
+3815pfp| - 43154p^p^ - 43154pjp| + 8192p| + 1582080pfp^ - 501760p?p| 
+65536pipi + 1582080plpi - 501760j3i;j3i + 65536j9lj32) 

TABLE 15: Polynomial expressions for (n) in the general case (2D periodic square planar lattices 
of increasing size) 
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N 




(n) 
Phc 


4 




.710835 


16 




.842153 


36 




.887094 


TABLE 16: 
N). 


Critical values p^^l (2D square planar lattices of increasing size with an even value of 


N 


(v) 




4 


(2/3)(2176j9ij52 - 2082pIpI + 923j5fj5| - 2082pfp| + UASp^pi 
+1448pfp2 - 320pf - 320p| + 1088pf - 4864pfp2 + 1088p| - 4864pip| 
+6418p?pi)/(13p?pi - 22pip| + 8pl - 22pfp2 + 16piP2 + 8plf 




9 


8(1312pij32 - 239p?pf + 77pIpI - 239pfp^ + 232p^pi + 232pfp2 
-72pl - 72p^ + 656pf - 1832pfp2 + 656p| - 1832pip^ + 1464pfp|)/ 
{-16pIp2 - 16piP2 + ^PipI + + 16piP2 + 8^2)^ 




16 


(8/135) (44708352^?^P2 + 44708352p^pi + 1 879900 160p?pi + 1409925120pfp| 





-7866368000pfpi + 1409925120pfpf - 7866368000p?pf + 16208510208pfpf - 3950903296pfpi 
-5062656pJ + 10880187392pfp^ - 14811384576pfpi + 10880187392pf - 14811384576pfpf 
+10174354048pfp^ + 563970048pfpi + 563970048pfp2 + 93995008^^ + 93995008p^ 
-5062656p^ - 3094999904pfpf + 5409620192pfp^ - 5210329472pfpf + 2820640640pfp^ 
+869041300p5p2 - 3094999904j3?p| + 5409620192p?p| - 5210329472p?p| + 2820640640p?j3| 
-804356096j9^j3i - 804356096pfp2 - 366246240j3fp^ + 23928810Gpf j?^ - 82905852pf j?^ 
-82905852p^p6 + ii685177p^p5 + 239288100p^p| - 366246240pjp^ + 322443168p^p? 
-3950903296pfpl + 322443168p^pJ - 163920768p]'p| - 163920768plpl) / {U7pfp^ 
-664pfp| + lOOOpfpf - eOSpjpi + 128pf - 664pfp^ + 2000pfpf - 1824pfp| 
+512p^pi + lOOOpfpl - 1824pIpI + 768plpl - 608pfp2 + hl2p\p2 + I28p\f 



TABLE 17: Polynomial expressions for {v) 


in the general case (2D periodic square planar lattices 


of increasing size) 




N 


{v) 
Pl'c 


4 


.702484 


16 


.837365 


36 


.885175 



TABLE 18: Critical values pfl (2D periodic square planar lattices of increasing size with an even 
number of sites N). 
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